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ABSTRACT 

A strong case has been made that radio waves from sources within about half 
a degree of the Galactic Center undergo extreme diffractive scattering. However, 
problems arise when standard ( "Kolmogorov" ) models of electron density fluctua- 
tions are employed to interpret the observations of scattering in conjunction with 
those of free-free radio emission. Specifically, the outer scale of a Kolmogorov 
spectrum of electron density fluctuations is constrained to be so small that it is 
difficult to identify an appropriate astronomical setting. Moreover, an unaccept- 
ably high turbulent heating rate results if the outer scale of the velocity field 
coincides with that of the density fluctuations. We propose an alternative model 
based on folded magnetic field structures that have been reported in numerical 
simulations of small-scale dynamos. Nearly isothermal density variations across 
thin current sheets suffice to account for the scattering. There is no problem 
of excess turbulent heating because the outer scale for the velocity fluctuations 
is much larger than the widths of the current sheets. We speculate that inter- 
stellar magnetic fields could possess geometries that reflect their origins: fields 
maintained by the galactic dynamo could have large correlation lengths, whereas 
those stirred by local energetic events might exhibit folded structures. 

Subject headings: Galaxy: center — ISM: general — radio continuum: ISM - 
scattering 



1. Introduction 

It appears that radio waves from sources within about half a degree of the Galactic 
Center (GC) undergo extreme diffractive scattering. Observations of SgrA* (Davies, Walsh 
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& Booth 1976; Rogers et al. 1994) and of maser spots in several OH/IR stars (van Langevelde 
et al. 1992; Frail et al. 1994) have established that the angular broadening is o b s ~ 1" at 
v ~ 1 GHz. If the scattering region is located more than a few kpc from the GC, this would 
correspond to enhanced, but not unusually large, levels of scattering; NGC 6334B is the 
record holder with a scatter broadened disk of size ~ 7" at v ~ 1 GHz (Moran et al. 1990). 
However, if the scattering region is located close to the GC (say, within 1 deg, corresponding 
to a distance R xs 150 pc from the GC), then the GC region would be a site of extreme 
scattering, van Langevelde et al. (1992) argue that, if the scattering region is also responsible 
for the free-free absorption toward the GC, then upper limits on the optical depth (Pedlar et 
al. 1989; Anantharamaiah et al. 1991) constrain R to be in the range 0.85 — 3 kpc. In more 
recent work, Lazio & Cordes (1998, hereafter LC98) make the case that R ~ 130 pc. Since 
the Sun's distance from the GC is ~ 8.5 kpc, this would imply that the angular scattering is 
larger than the observed angular broadening by a factor 8500/130 ~ 65: 

fl scat w 65 TZ- 1 9 ohs ~ 3 x 10- 4 ft-V 9 - 2 , (1) 

where we have set R = 130 TZ pc, and u 9 = z//(GHz). Spergel & Blitz (1992) estimate the 
gas pressure in the Galactic bulge to be p = nT ~5x 10 6 Kcm~ 3 . Studies based on X-ray 
emission from the GC (Muno et al. 2004) suggest that, even on the smaller spatial scales 
of interest to us, the gas pressure may not be very different. A pressure of 5 x 10 6 Kcm~ 3 
is about 10 2 ' 5 times higher than in the local interstellar medium, hence it may not seem 
unreasonable to expect the GC to be a region of extreme scattering for radio waves. However, 
as LC98 demonstrate, problems arise when standard models of electron density fluctuations 
are employed to interpret the observations of scattering in conjunction with those of free-free 
radio emission. 1 



2. Problems With The Standard Interpretation 

LC98 estimated the brightness temperature of free-free emission toward five highly scat- 
tered OH maser sources near the GC, based on the 10 GHz survey of Handa et al. (1987). 
From Table 3 of LC98, we derive a mean brightness temperature, ~ 0.25 K at v — 10 GHz . 
Since the region is optically thin, the optical depth tr ~ (T&/T) ~ 2.5 x 10 _5 T 4 _1 , where 
T 4 = T/(10 4 K) is the gas temperature. Expressed in terms of the emission measure (EM), 
Tff 5 x 10~ 7 i>g 2 T 4 3 ^ 2 EMi, where EM\ = _EM/(pccm~ 6 ), and we have set the averaged 



1 Similar, but less severe, problems arise for NGC 6334B which lies behind a galactic HII region (Moran 
et al. 1990). 
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Gaunt factor g s p=* 10. Hence, 

EM ~ 5 x 10 3 T 4 1/2 pc cm" 6 . (2) 

The emission measure contributed by a medium with pressure p = 5 x 10 6 JFK cm -3 , is 
EM ~ /n 2 i? ps 3.3 x 10 7 / T 2 T 4 ~ 2 K pccm" 6 , where / is the volume filling factor. Provided 
all the gas that contributes to the EM is at temperature T4, / ~ 1.5 x 10~ 4 T 4 ^ 2 JF -2 Tlr x 
and the mass of this gas is M ~ 1.5 x 10 4 T 3/2 T' x ll 2 M Q . Obviously the constraint / < 1 
implies T 4 < 30 T^K 2 ^. 

In the standard interpretation of radio wave scattering in the interstellar medium, the 
electron density fluctuations are assumed to follow a Kolmogorov spectrum, (Ane/n) ps 
(£/L) 1/3 for £ mhl < £ < L. For SgrA*, # scat oc v~ 2 for v 9 < 30. The V 2 scaling implies 
that the scales contributing dominantly to strong, diffractive scattering are smaller than 
the inner scale: £ ni \ n > A/# scat , where A = cjv is the wavelength of the radio waves. 2 
Because \/0 scat oc v>, we obtain a lower limit on £ min which is higher at higher frequency. 
Using equation (1), with v g = 30, gives £ min > 3 x 10 6 7£cm. The scattering angle, 6> scat ~ 
(z/ p /z/) 2 (/i?/L) 1 / 2 (L/£ min ) 1 / 6 , where v v m 10 4 (n/ cm -3 ) 1 / 2 Hz is the plasma frequency. In this 
expression for ^ scat , we can write fR ^ EM/n 2 . Using equations (1) and (2) for # scat and 
EM, the condition £ min > 3 x 10 6 K cm implies that, 

L<3x 10" 8 T 3/4 ^ 5/2 pc. (3) 

The upper limit to L is remarkably small and independent of the pressure. 3 Moreover, if 
only a fraction F of the plasma that contributes to the free-free emission is responsible for 
the scattering, L would be smaller by a factor F 3 ^ 2 . This, in essence, is LC98's argument. 
Henceforth we set T — 1Z — 1. 

The choice of a small outer scale, constrained by equation (3), introduces two problems. 
As was realized by LC98, it implies an unacceptably high heating rate if L is also the 
outer scale of the velocity fluctuations. Dissipation by Kolmogorov turbulence would heat 
the gas at rate, t^ rh ~ c s /L > 10~ 5 T 4 1 ^ 4 s _1 , whereas radiative cooling would be much 
slower, t~* ol ~ 3 x 10 -7 T / f 2 A_ 2 i s _1 (here c s ~ lOT^kms -1 is the sound speed, and A_ 2 i = 



2 If scattering was dominated by scales between £ m ; n and L, we would have # sca t oc v 2 2 . 

3 The problem of a small outer scale for density fluctuations arises more generally. Anantharamaiah & 
Narayan (1988) estimated that enhanced scattering in the inner galaxy required fractional density fluctua- 
tions ~ 10 on the outer scale. For a Kolmogorov spectrum, the same enhancement may be achieved with 
order unity fractional density fluctuations, with an outer scale that is 10 3 times smaller. In this case, the 
outer scale of the velocity fluctuations can be much larger than the outer scale of the density fluctuations. 
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A/(10~ 21 ergs -1 cm 3 ); A_ 2 i ~ 1 when T 4 ps 10). It is also difficult to identify an appropriate 
astronomical setting in which nonlinear density fluctuations on a scale as small as L might 
arise. We have considered radiative shocks and interfaces between neutral and ionized gas 
formed by ionizing radiation or hot gas incident upon the surfaces of a molecular clouds. 
Even at the high pressure in the GC region, each case yields a length scale that is at least a 
few orders of magnitude larger than L. 



3. Folded Fields As A Solution 

Folded magnetic field structures appear in numerical simulations of small-scale, turbu- 
lent dynamos (Schekochihin et al. 2004). Of relevance here are the electrical conductivity, 
a ~ 10 13 T 4 3 ^ 2 s -1 , and the kinematic viscosity, z/ v i S ps 2 x 10 15 T 4 7 ^ 2 cm 2 s" 1 , for plasma with gas 
pressure, p ~ 5 x 10 6 K cm -3 . 4 The magnetic Prandtl number, Pr m = 47rc -2 o"z/ v i S Ri 3x 10 8 Tl, 
is large so dissipation of velocity fields occurs on much larger scales than dissipation of mag- 
netic fields. Schekochihin et al. (2004) propose a model for the organization of magnetic fields 
that appears attractive; we describe this briefly below, before considering its implications 
for density fluctuations. 

Let an incompressible fluid, permeated by a weak magnetic field, be stirred on an 
outer scale L, with random velocity Vl- Within a few stirring times, r L ~ (L/vl), the 
kinetic energy cascades turbulently, and creates velocity fluctuations on smaller spatial scales 
through nonlinear hydrodynamic interactions. The rms velocity across a scale i may be 
expected to follow a Kolmogorov spectrum, vg ~ vl^/L) 1 ^ 3 , for £ vis < £ < L, where the 
inner scale is £ V - 1S ~ L(u vis / Lvl) 3 ^ 4 <C L. Eddies of scale i turn over on times, ~ (£/ve) ~ 
tl(£ I L) 2 / 3 . Hence the early evolution of the weak magnetic field will be dominated by the 
stretching action of the smallest eddies, of size ~ £ v ; s , because their turn-over time is the 
shortest. These eddies cease to be effective at stretching the field lines when the magnetic 
energy density becomes comparable to their kinetic energy density, i.e. when B 2 ~ m p nv 2 is . 
Then larger and more energetic eddies take over and the magnetic energy density continues 
to grow until it achieves approximate equipartition with the kinetic energy of the largest 
eddies: B 2 ~ m p nv 2 L . It is the geometry of the magnetic field that is of particular interest to 
us. According to Schekochihin et al. (2004), it has a folded structure, with parallel correlation 
length ~ L: after a distance L, a typical field line reverses direction sharply, and folds back 
on itself. Sheets of direction-reversing folded fields of thickness d are separated by current 



4 These scalar transport coefficients only apply in directions parallel to the local magnetic field. 
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/ 2 \ 1/2 

/ C \ ^10^.-3/4 1/2 



sheets of similar thickness, 5 

^U^ TL J ~ 10lUT ^ /V 6 /Zcm ' W 

where Tq = tl/(1 Myr) . For an isothermal gas at fixed total pressure, the fractional density 
perturbation across a current sheet of thickness d is, 

*2 „ 3 (5) 

n SnnkT f V ; 

Next we estimate scattering of radio waves by folded fields. Consider the idealized case 
of an ensemble of plane-parallel current sheets, each of radius L and thickness d <C L. 
Let the sheets fill space statistically homogeneously with filling factor, /, and be oriented 
randomly. 6 The rms phase difference across a transverse scale d, accrued along a path length 
R ^> L, can be estimated by imagining a ray going through (R/L) independently oriented 
plates, of which only a fraction (d/L) is oriented favorably enough to each contribute to the 
phase a path length ~ L: 

For a more realistic case, we take the sheets to have radii of curvature r c , perhaps caused by 
Alfven waves which can propagate along folded fields. To estimate the rms phase difference 
for r c > L, there are two limits to consider. For L < \Jdr c , the situation is similar to that 
for plane-parallel sheets, and A$ is given by equation (6). For \fd~r~ c < L < r c , a ray goes 
through {R/L) independently oriented plates, of which a larger fraction {L/r c ) is oriented 
favorably enough for each plate to contribute path length ~ \fd~F c . Therefore, the expression 
for A$ remains unchanged, so long as r c > L. We can also arrive at the above conclusions 
more formally by estimating the phase structure function, due to isotropically oriented sheets 
(see Appendix). The angular scattering, 6> scat ~ (A/d)A$ ~ {v p /v) 2 {fR/d) l / 2 f3~ l ) depends 
on L through d but is otherwise independent of L. We write fR w (EM/n 2 ), and use 
equation (4) for d, to obtain a general expression for the angular scattering: 

scat - (7) 2 (^) V2 p- 1 - 2 >< 10 ~ 6 17 Vr 1 • (7) 



5 Foldcd fields in cartoon form, and as produced in a simulation, are shown in Figs. 10 & 15 of Schekochihin 
ct al. (2004). 

6 A modest non-random orientation could account for the anisotropic images of scatter broadened sources 
near the GC. 
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This is our main result. It is independent of the geometrical distribution of the scattering 
material, but this deserves further attention. 

The scattering material must cover most of the area as seen from outside R. Let us 
suppose that folded fields fill an approximately spherical shell of radius R and thickness 

AR ~ fR ~ 2 x 1(T 2 T 4 5/2 T~ 2 pc . (8) 

Assuming that turbulent stirring occurs at sonic speeds, vl ~ lOT^kms -1 , we find 

r 6 ~ 2 x 10- 3 T| (-^-) . (9) 



Substituting this value of r 6 and equation (2) for EM in equation (7), we obtain, 

«-~ 5 x 10 (wj {-) ' (10) 

which compares well with the value of # sca t, derived from observations, given by equation (1). 
The thickness of the current sheets can be estimated by substituting equation (9) for r 6 in 
equation (4): 

, / T \ 1/2 

d~ 6 x 10 8 Tj/ 4 cm. (11) 

The scattering is strong because 

A»~£u~l*(g)(A) V \ (12) 

One motivation for considering scattering from a shell of warm ionized gas is the presence 
of a lobe of emission surrounding the GC, discovered by Sofue & Handa (1984) in a 10 GHz 
survey. Comparison with the 5 GHz survey of Altenhoff et al. (1978) enabled them to decom- 
pose the emission into thermal and nonthermal components, and establish that the thermal 
component arises mostly from a shell-like feature (Sofue 1985). This is one of the larger of 
the many sources in the GC region (LaRosa et al. 2000), with an angular size exceeding a 
degree. Even so, its radius ~ 80 pc appears to be only a little more than half the value of 
~ 130 pc we assumed, based on LC98's location of the scattering region. This implies that 
AR might be only a little larger than half the value given in equation (8). However, our 
estimates of 6* scat and A$ given above are hardly affected, because of their weak dependence 
on AR. 
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4. Summary 

Turbulence giving rise to folded fields may account for extreme scattering in the GC 
region because large fractional density fluctuations occur across the thickness of the current 
sheets, d, which is much smaller than the outer scale, L. Although the current sheets are 
thin, over-heating is not a problem. The ohmic dissipation rate per volume due to current 
sheets is ~ (cB/d) 2 a L ~ B 2 /t l . If, as we have argued, the magnetic energy achieves near- 
equipartition with the kinetic energy, then B 2 ~ m p nv 2 Ll and the ohmic dissipation rate is no 
larger than the turn-over rate of the largest eddies: t~^ m ~ t^ 1 , which is much smaller than 
t~* ol . The large cooling rate also ensures that the plasma behaves nearly isothermally, an 
assumption used in our estimate of density fluctuations across the thickness of the current 
sheets, given in equation (5). 

In our discussion of radio-wave scattering we assumed that the current sheets are ori- 
ented isotropically. It then proves convenient to formulate radio-wave scattering in terms 
of angle-averaged quantities. The effective, isotropic power spectrum of electron density 
fluctuations turns out to be shallow, being proportional to k~ 2 (see equation A2). Such a 
spectrum gives rise to large density fluctuations on spatial scales ~ d, but contributes little 
to density fluctuations on large spatial scales. This difference provides an observational test 
of the folded field hypothesis. Where diffractive and refractive scintillations have been de- 
tected in the same source, their relative magnitudes suggest that a Kolmogorov spectrum of 
density fluctuations extends up to at least the refractive scales. However, refractive scintilla- 
tions are generally not detected in strongly scattered sources. This does not conflict with the 
Kolmogorov model which predicts them to be both weak and slow, but it is also compatible 
with the folded field hypothesis. Higher frequency observations could provide a decisive test. 
For a Kolmogorov spectrum, the fractional flux modulation, m, is of order the cube-root of 
the ratio of the diffractive scale to the Fresnel scale. With parameters appropriate to SgrA*, 
we find m ~ lO^TZ^u^ 2 . 

Magnetic fields might possess different geometries in different regions of the interstellar 
medium, depending on their origins. Some regions could possess strong, mean magnetic 
fields, correlated over large distances, generated, perhaps, by large-scale dynamos. These 
could be the sites of anisotropic Kolmogorov turbulence (Goldreich & Sridhar 1995), re- 
sponsible for the general level of diffractive and refractive interstellar scintillation. There 
could be other sites in the interstellar medium, permeated with folded fields, generated by 
small-scale dynamos, contributing to extreme diffractive scattering, but little to refractive 
scattering: the Galactic Center could be one such region. 

We would like to inject a cautionary concluding note. Our model of scattering applies 
results on the growth of magnetic fields in small-scale dynamos from Schekochihin et al. 
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(2004). These are based on MHD with scalar diffusivities. It is an open question as to 
whether they can be applied to the low density, magnetized plasma in the GC region. We 
expect that our estimate for magnetic diffusivity, which sets the thickness of the current 
sheets, is valid. But the reduction of the kinematic viscosity in directions perpendicular to 
the magnetic field is a concern. It is possible that folded fields would unwind so rapidly that 
they could not be maintained. 



A. Diffractive scattering by an ensemble of current sheets 

Here we offer a physical derivation of equation (6), by estimating the phase structure 
function due to randomly distributed and isotropically oriented current sheets. A useful 
intermediate step is to calculate C(r), the effective isotropic density-density correlation 
function on separation r. Its Fourier transform is P(k), the effective isotropic power spectrum 
of density fluctuations. Consider an emsemble of flat current sheets, each of thickness d 
and radius L ^> d, all oriented in the same direction with unit normal n. The density- 
density correlation function is significantly non-zero only when the separation, r, is such 
that |r ■ n| < d , and r± = |r — n(r • n) | < L . When averaged over all directions of n, we will 
obtain C(r). Equivalently, we may keep n fixed, and average over all directions of r. When 
d < r < L, the solid angle of the region of intersection between the plate-like region (of 
radius L and thickness d), and a concentric sphere of radius r, is ~ (dr/r 2 ) ~ (d/r). Hence, 

C(r) ~ (An) 2 ^; for d<r<L. (Al) 
The corresponding effective isotropic power spectrum is "shallow" : 

P(k) ~ (An) 2 ^- ; for L" 1 < k < d~\ (A2) 
k 

P(k) ~ (An) 2 L 2 d when kL < 1 and P(k) ~ for kd > 1. Let us consider the case when the 
current sheets are curved, with radius of curvature, r c > L. The solid angle of the region of 
intersection between the curved current sheet-like region and a concentric sphere of radius r 
(in the space of separations) is, ~ (dr^/r 2 ). We estimate, r± ~ r (1 — r 2 /8r 2 ) ~ r, because 
r/r c < L/r c < 1 . Therefore the solid angle of intersection is (d/r), and the expressions for 
C(r) and P{k) are as given above. As may be seen, the dominant contribution to density 
fluctuations comes from scales close to d. 

In the thin-screen model of radio-wave scattering, the entire effect of the interstellar 
medium is specified by the (gaussian) random phase pattern imprinted on a wave front, 
as it passes through a "phase screen", placed between the source and the observer. The 
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statistical properties of the screen are completely described by the phase structure function, 
-D(s), which is defined as the mean square phase difference across transverse separation s on 
the screen. When the power spectrum is isotropic, D depends only on s = |s|: 



" Jo 

where r e = e 2 /m e c 2 is the classical electron radius. For the shallow spectrum of equa- 
tion (A2), it is straightforward to make the estimate, 



where A$ is given by equation (6). If A$ < 1 , the scattering is weak, whereas a range of 
scales (smaller than, and of order d) can contribute to strong scattering, when A$ > 1 7 
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